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Abstract. The Medusa (NGC 4194) is a well-studied nearby galaxy with the disturbed appearance of a merger and evidence 
for ongoing star formation. In order to test whether it could be the result of an interaction between a gas-rich disk-like galaxy 
and a larger elliptical, we have carried out optical and radio observations of the stars and the gas in the Medusa, and performed 
W-body numerical simulations of the evolution of such a system. We used the Nordic Optical Telescope to obtain a deep V-band 
image and the Westerbork Radio Synthesis Telescope to map the large-scale distribution and kinematics of atomic hydrogen. 
A single Hi tail was found to the South of the Medusa with a projected length of ~ 56 kpc (~ 5') and a gas mass of 7 • 10** Mq, 
thus harbouring about one third of the total Hi mass of the system. Hi was also detected in absorption toward the continuum in 
the center. Hi was detected in a small nearby galaxy to the North-West of the Medusa at a projected distance of 91 kpc. It is, 
however, unlikely that this galaxy has had a significant influence on the evolution of the Medusa. The simulations of the slightly 
prograde iirfall of a gas-rich disk galaxy on an larger, four time more massive elliptical (spherical) galaxy reproduce most of the 
observed features of the Medusa.Thus, the Medusa is an ideal object to study the merger-induced star formation contribution 
from the small galaxy of a minor merger. 
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1. Introduction 

The Medusa galaxy (NCG4194, Arp 160, IZw33) is a nearby 
galaxy (D = 39 Mpc) with the disturbed appearance typical 
of a merger (see Fig.[T]i. It has been interpreted as a candidate 
for a minor merger between a large elUptical and a small spiral 
galaxy, a so-called S+E merger (Aalto & Hiittemeister"2000). 
The galaxy earned its nickname from the fuzzy and knotty op- 
tical tail to the North, reminiscent of hair or tentacles. A sharp 
shell-like structure is visible on the opposite side of the galaxy, 
confining the optical main body. A second, fainter shell can be 
seen further out on deep images (e.g., Jtitte et al.|2008| l. 

Major mergers between two gas-rich disk galaxies can lead 
to ultraluminous infrared galaxies (ULIRGs, L/r/s > 10'^ Lq) 
due to an extreme increase of star formation ( [Sanders et al.| 
2003 I. Strong gas inflows are triggered in these mergers, in- 
creasing the gas density in the centre and thus leading to the ob- 
served starburst (e.g., [Barnes & Hernquist|1991|[l996||Mihos| 
& Hernquist 19961. In these simulations the merger remnant 



evolves into a regular elliptical galaxy. In contrast, in minor, 
i.e., unequal-mass mergers only a moderate, if at all, starburst 
is induced implying a lower far-infrared (FIR) luminosity (less 
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than lO'^ Lq) ( Cox et aT]|2008 1. The triggered star formation 
and the remnant galaxy type strongly depend on the mass ra- 
tio of the merging pair, but not so much on the total gas mass 
Poumaud et al.[|2005] [Cox et al.|[2008| l. In contrast to most 
nearby minor mergers, in which the major galaxy is a spiral and 
thus dominating the star formation (e.g., M51), the Medusa 
gives us an ideal opportunity to investigate the influence of 
the merger on the smaller companion because it is the only 
gas-rich paitner here. Since the merger-induced star formation 



does not depend on the gas content of the companion (Martig 
& Bournaud|2008 1, the results of the Medusa might be appl 



cable to other minor mergers, independent on the type of the 
large partner 

Another significant difference between ULIRGs and galaxies 
like the Medusa is the extent of the starburst induced by the 
merger While the burst region in ULIRGs is compact and typ- 
ically concentrated to the inner kiloparsec, (e.g 



Sanders & 



Mirabel] 1996 1, the Medusa exhibits an extended region of on- 



going star formation which is, however, not as intense as in 
ULIRGs (|Armus et al.|1990| [Prestwich et al.|1994]i.|Weistrop[ 



et al. (2004 1 analyzed optical and ultraviolet Hubble Space 



Telescope (HST) images of the nucleus and found star-forming 
knots younger than 20 Myr in the center They derived a star 
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formation rate (SFR) of ~ 6 M© yr"', but argued that the over- 
all SFR in the Medusa might be as high as 30 Mg yr"' since the 
knots produce only 20% of the UV flux in the center. Based on 
recent HST spectroscopic observations, Hancock et al. ( 2006| l 



Table 1. Basic properties of NGC4194. The distance is based 
on Ho - 75 km s"' Mpc"'. 



estimate an even higher total SFR of ~ 46 Mq yr^ and ages of 
individual star forming regions of 5.5-10.5 Myr 

The Medusa possesses a reservoir of gas, both molecu- 
lar and atomic. High resolution aperture synthesis maps ob- 
tained with the Owens Valley Radio Observatory revealed ex- 
tended CO emission on a total scale of 25" (4.7 kpc) ( |Aalto 
|& Hutt emeister|2000p . The CO morphology is complex, occu- 
pying mainly the center and the north-eastern part of the main 
optical body. The extended CO gas follows two prominent dust 
lanes: one which is crossing the central region at right angle 
with respect to the optical major axis, and a second which 
curves to the north-east and into the beginning of the north- 
ern optical tail. |Aalto & Hiittemeister ( 2000| l suggested that the 
central starburst is being fueled by gas flows along the central 
dust lane. The presence of multiple spectral peaks (AV « 200 
km s~') suggests that at least two velocity systems coexist in 
the CO emitting gas. Either this is an effect of substantial warp- 
ing of the main body, or simply the unrelaxed, overlapping sys- 
tems of the progenitor galaxies. The multiple features seem to 
be related to the extended dust lanes, but several velocity com- 
ponents co-exist also at the base of the optical tail. Single-dish 
measurements revealed the presence CO emission even inside 
the optical tail out to 4.7 kpc away from the center ( Aalto et al. 
2001)1 whe re there is no sign of ongoing star formation ( Armus 



et al.|1990 1. The total molecular gas mass in the inner 2 kpc was 
estimated to 2- 10^ Mq. Thus, the star formation efliciency is al- 
most as high as in ULIRGs although the molecular gas density 



is significantly lower. (Aalto & Hiittemeister 2000). Besides 
uncertainties in the molecular gas conversion factor, this gives 
a hint that the ISM in these mergers is difi'erent compared to 
the molecular cloud complexes in normal galaxies and also in 



extreme starbursts. | Aalto & Hiittemeister (20001 argue that the 
observed solid-body rotation may off'er a shear-free environ- 
ment. This may support cloud collapse and thus star formation. 
The star formation rate might have increased because a large 
amount of molecular gas was transported to the centre on short 
time scales by the interaction. 

The galaxy contains a similar amount of atomic gas as 
molecular gas {Mhi~ 2.2 • 10^ Mq), inferred from the single- 



dish Hi emission spectrum obtained by Thuan & Martin ( 1981 1 
using the Greenbank 91m telescope. Atomic gas has recently 
been detected in absorption toward the center of the Medusa, 



through sub-arcsecond MERLIN observations ( Beswick et al. 



|2005| l. The kinematics of the Hi gas in the central region of 
the galaxy is similar to that of the molecular gas traced by the 



CO emission ( jAalto & Hiittemeister 2000 1, which suggests that 
the atomic and molecular component of the cold interstellar 
medium are probably physically related. Two compact radio 
sources were found in the central parsec, embedded in a region 
of difiiise radio continuum emission. Even though the presence 
of a weak AGN cannot be ruled out completely, most (if not all) 
of the continuum emission is related to the starburst ( [Beswick] 
let al.|2005| l. 



RA (2000) 

DEC (2000) 

Vhei.opt (kms"') 

D (Mpc) 

LedO'Lo) 

Lf,r(10'"Lo) 

FIRfluxat60yum(Jy) 

FIR flux at lOOyum (Jy) 



12 14 09.5 
+54 31 37 
2442 
39 
15.7 
8.5 
25.66 
26.21 



Table 2. Observing parameters of the WSRT observations. 



observing parameters 
date 

central frequency 
total bandwidth 
number of channels 
velocity resolution 
synthesised beamwidth 
primary beam 
primary calibrator 



12.7.2003 
1408.6 MHz 
20 MHz 
256 
4.12 km S-' 
23" X 19" 
36' 
3C147 



NGC 4194 displays conflicting evidence about its merger 
age. The sub-arcsecond resolution MERLIN data revealed two 
compact radio components in the center, separated by ~ .35, 
corresponding to 65 pc only ^Beswick et al.|[2005 l. The opti- 
cal tail is fairly diffuse, which is suggestive of an advanced 
merger where the progenitors have coalesced. In contrast, the 
fairly widespread CO distribution, irregular optical isophotes, 
multiple spectral CO features, and the distorted central body 
hint toward a younger system. 

In this paper, we present observations and A^-body nu- 
merical simulations of the stars and the gas in the Medusa 
merger We mapped the large-scale distribution and kinematics 
of the Hi gas with the Westerbork Radio Synthesis Telescope 
(WSRT) and used the Nordic Optical Telescope (NOT) to ob- 
tain a deep V-band image of the distribution of the stars in 
the Medusa. Then we carried out numerical simulations which 
we compared to the observations. Hi 21 cm line observations 
are a powerful tool to trace the dynamical history of interact- 
ing and merging galaxies and provide constraints to dynamical 
models of the history of the interaction. They complement CO 
millimeter-wave fine observations that trace mainly the inner 
part of galaxies where the molecular gas resides. Furthermore, 
they offer is the only way to trace the distribution and kine- 
matics of gas in large scale tidal features. The observations are 
presented in section 2 and 3, and the simulations in section 4. 
Follows a discussion. 

2. Observations and data reduction 

2.1. Hi observations 

The observations were caiTied out with the WSRT in July 2003. 
We observed one full 12h track in the maxi-short array config- 
uration under good conditions. As a flux calibrator, the stan- 
dard source 3C147 was observed before and after the target. 
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Fig. 1. Left: Nordic Optical Telescope V-band image of the Medusa. Right: Contour map of the 20 cm continuum emission 
superimposed on the V-band image. The synthesised beam is shown in the lower left corner. The contour levels are 0.2, 0.5, 2, 5, 
20, 50, 70mJy beam-'. 



No phase calibrator was observed, since the standard procedure 
at WSRT for phase calibration is to obtain self-calibration. In 
Tab.|2] we list the observing parameters and give information 
about the setup used. We chose parameters which allowed us 
to both cover a large velocity range, because Hi lines can typi- 
cally be broader in mergers than in normal disk galaxies, and to 
reach a high resolution in velocity to be able to detect narrow 
features. 

The data were reduced using MIRIAD together with addi- 
tional, WSRT specific routines written by T. OosterlocQ The 
quality of the data was good, so that only little flagging of in- 
terferences was necessary. First, the data had to be calibrated 
to system temperature scale. Since no secondary calibrator was 
observed, after flux and bandpass calibration using the primary 
calibrator, a continuum image was created to perform self- 
calibration on continuum sources within the field. To achieve a 
deep CLEANing, we iteratively created a mask which includes 
emission only and which defines the regions to be CLEANed. 
After CLEANing and self-calibration a new mask was created 
and the process repeated until no improvement could be made. 
The results were then applied to the target galaxies. 

CLEANing the Hi cube was done in a similar way to 
CLEANing the continuum. After continuum subtraction, we 
applied Hanning smoothing, since this was not done during the 
observations. Again, we used as mask containing emission only 
to define CLEANing boxes. This procedure was repeated 3 
times until the noise level was reached. Integrated intensity and 
velocity maps were created using the MOMENT task in MIRIAD. 
Finally, primary beam correction was applied to all data. 

We produced two sets of maps: naturally weighted ones 
(with an rms noise level of 0.4 mJy beam resolution of 
31" X 27"), and maps obtained using a ROBUST factor of 0.5, 

' http://www.astron.nl/~oosterlo/wsrtMiriad/ 



which gives a better resolution with slightly less sensitivity 
(noise level of 0.6 mJy beam"', resolution of 23"x 19"). Since 
the morphological features were the same in both sets of maps, 
we present here only the ROBUST . 5 maps in which structures 
in the Hi tail appear more clearly. 

2.2. Optical imaging 

We also present here optical Johnson V-band images taken with 
ALFOSC mounted at the Nordic Optical Telescope (NOT) on 
La Palma, Canary Islands in January 2003. The field of view 
of ALFOSC is 6'A x 6f4, thus large enough to cover not only 
NGC4194 completely but also a sufficient region of sky for 
a proper background subtraction. The integration time was 10 
minutes. Standard reduction was applied to the data using the 
IRAF software. A further description of the observations and 
reduction as well as a detailed analysis of the optical data will 
be given in Jiitte et al. ( 2008[ l. 



3. Observational results 

3.1. Radio continuum 

Figure[T]shows the NOT V-band image of NGC 4194 on which 
contours of the 20 cm radio continuum emission have been su- 
perimposed. The flux density of the source is 0.108 Jy. We see 
continuum emission extending out to a deconvolved size of 6'.6 
(1.2kpc) associated with the galactic main body. This is not 
surprising, since the 20 cm continuum flux is a tracer of star 
formation and NGC 4194 exhibits an extended region of in- 
tense star formation. The continuum map shows a smooth dis- 
tribution center at the position of the optical nucleus. 
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Fig. 2. The Hi distribution as contours overlayed on an optical image of the field around NGC4194 from the Digitized Sky 
Survey. The contour levels are 0.01,0.05,0.1,0.15,0.2,0.3,0.4,0.5 Jy beam km s"'. The bar in the lower right corner marks a 
projected distance of 10 kpc, assuming a distance of 39Mpc. The lack of contours at the center of NGC4194 is due to an 
absorption source. The beam is shown in the lower left corner. 



3.2. Neutral hydrogen 

Figure|2]shows the integrated Hi intensity map overlayed on an 
optical image from the Digitized Sky Survey (DSS). Hi gas is 
found all over the optical extent of the Medusa. Remarkable is 
the single long Hi tail to the South. Hi is also found in absorp- 
tion toward the center. We also detect Hi associated with the 
small nearby galaxy SDSS J121326.03+543631.7 identified in 
the Sloan Digitized Sky Survey about 8' to the north-west. 

The Hi properties of NGC4194 and of the nearby dwarf 
galaxy are summarised in Tab. |3] The integrated Hi spectra are 



shown in Fig. |3] We calculate a total Hi mass of NGC4194 
of 2 ■ 10^ Mq, in good agreemen t with the Hi mass determi- 
nation by |Thuan & Martin| (|l98l') of 2.2 ■ 10'^ Mq. This shows 
that our interferometric observations do not miss any extended 
emission. 



3.2.1 . The main body of NGC 41 94 

Figure|4]shows the Hi distribution, the velocity field and the ve- 
locity dispersion map. The channel maps are shown in Fig.|5] 
Figure|6]shows the Hi integrated emission as isocontours over- 
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Fig. 3. Integrated Hi spectra of NGC 4194 (left) and the companion dwarf galaxy to the North- West (right). 



Table 3. Hi properties of the Medusa merger and of the nearby 
dwarf galaxy SDSS J121326.03+543631.7. The values listed 
for NGC 4194 include the contribution from the Hi tail. 



Hi properties 


NGC 4 194 


Tail 


Nearby dwarf 


Vccnt (km/s) 


2500 




2370 


range (km s"') 


2400 - 2580 


2470 - 2520 


2350 - 2390 


extension (kpc) 


81 


56 


14 


Fm (Jykms"') 


5.65 


2.03 


0.20 


M^KlO'Mo) 


2.0 


0.7 


0.07 


MhiI Lb 


0.18 







laid on our optical V-band image, which is deeper than the DSS 
image but covers only the upper part of the Hi tail. The Hi 
emission peaks north of the absorption feature, aligned with 
the beginning of the optical tidal tail and some bright knots, 
which might be proto-star clusters. The Hi embeds the optical 
tail completely, but is not considerably extended further out, 
i.e., no pronounced Hi tail as a counterpart of the optical tail is 
found. The optical tail seems to be embedded in a smooth Hi 
distribution somewhat extended towards the north-west, maybe 
coincidentally in the direction where the possible companion is 
found. 

In the optical, we see two shells in the south, a bright one 
close to the center and a second, fainter one at a larger radius. 
We do not see any feature, neither in the Hi distribution nor in 
the velocity and velocity dispersion maps which is related to 
these shells. 

The velocity map (Fig. 4, bottom left) shows a slightly dis- 
turbed rotating disk, represented by a typical spider diagram, 
with regions of higher velocity on the eastern side and lower 
velocity on the western side on the main body of the Medusa. 

The Hi velocity dispersion is larger in the galactic body 
than in the tail. The highest values are measured north-east 
of the absorption region, where there is a slight dip in the 
Hi distribution. This region is also roughly related to the 



ity components are found, mimicing high velocity dispersion in 
an integrated dispersion map. It seems therefore likely that also 
the high dispersion in Hi is related to the dust lane rather than 
the starburst region, even though the spatial resolution of the 
Hi presented here is not sufficient for a detailed investigation. 

3.2.2. Hi in absorption toward the center 

Hi is seen in absorption against the underlying continuum 
source of the star forming nucleus. Fig. |7] shows the mea- 
sured spectrum. The multicomponent velocity structure clearly 
seen in the MERLIN data and discussed by |Beswick et al. 
( 2005[ l is tentatively seen here, too. Our spectrum shows at 



least two components, 
~ 2570 km s-i. 



one at 



2500 km s and one at 



dust lane where extended molecular gas was found (Aalto & 
Hiittemeister|2000 1. In the molecular gas maps multiple veloc- 



3.2.3. The Hi tail 

We detected a single remarkable Hi tail going to the South out 
to ~56 kpc from the center The Hi tail contains several clumps 
embedded in smooth, diffuse Hi emission. Neither in our V- 
band image (which, unfortunately, covers only the upper part 
of the Hi tail) nor in the larger DSS image did we find any op- 
tical counterpart of the Hi clumps. We determined the Hi mass 
in the tail as ~ 7 ■ 10^ Mg, i.e. 35% of the total Hi mass. The 
velocity in the tail is smoothly increasing, with the highest ve- 
locity of ~ 2510 km s"' at the tip of the tail. In the velocity 
dispersion map the smooth behaviour of the tail is found as a 
continuous but not disturbed increase of dispersion from the tip 
of the tail towards the galactic body. The dispersion is gener- 
ally higher at positions of the clumps, especially in two large 
clumps close to the galaxy's main body. 
The lowest velocity (~ 2455 - 2480 km s"') is found in the re- 
gion between the central Hi concentration and these two bright 
clumps, marking the beginning of the tail. In this region, there 
is a gradient in velocity from north-east to south-west, end- 
ing in an area having the lowest velocity in the system south- 
wards of the center. North-west of the two clumps the smooth 
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Fig. 6. Contours of the Hi intensity map overlayed on the V- 
band image. The optical image covers only the upper part of 
the Hi tail. The contour levels ai-e 0.01, 0.05, 0.1, 0.15, 0.2, 
0.3,0.4, 0.5 Jy beam-' kms"'. 

decrease of velocity is intercepted by a 'ridge' of slightly in- 
creased velocity. In the velocity dispersion map we see en- 
hanced dispersion there (~ 35 km s ' compared to ~ 10km s ' 
left and right of this 'ridge'), whereas in the Hi distribution 
map, if at all, there seems to be a slight dip in Hi intensity 
associated with this 'ridge'. 

3.2.4. Hi in a nearby dwarf galaxy 

We were also able to detect Hi in the dwarf galaxy 
J121326.03-H543631.7 at a distance of 8.06' (91kpc) to the 
North- West (a = 12''13"26.38% 5 = H-54°36'31"). Until now, 
it was not known that NGC4194 has a companion: our Hi 
measurements show that this dwarf has a central velocity of 
2380 km s"', which is only ~ 120km s"' lower than that of 
NGC4194, thus probably indicating that the two objects are 
physically connected. The measured Hi velocity is in good 
agreement with its optical velocity of 2436 km s"' (SDSS data 
release 2, 2004). We carefully checked our Hi channel maps 
for an Hi bridge, but even in the naturally weighted maps no 
connection is visible. Figure |3] shows the integrated Hi spec- 
tra of NGC4194 and the nearby galaxy derived from the data 
cube. With an Hi line width of ~ 60km s ' and an Hi mass 
of nnl ■ 10^ Mq this dwarf galaxy is significantly smaller than 
NGC 4194 and thus does not play an important role for the evo- 
lution ofNGC4194. 

4. N-body simulations 

We have performed numerical simulations in order to test 
whether the current appearance of the Medusa could be due 
to a gravitational encouter between a gas-rich disk-like galaxy 
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Fig. 7. The Hi absorption spectrum toward the center of 
NGC 4194. 



and an elliptical. Our goal was not to reproduce every detail 
of the Medusa, but rather to see whether a simple type of en- 
counter could generate the observed overall distribution of the 
stars and the gas, and in particular the following striking pat- 
terns: i) the diffuse region of stellar emission to the north (the 
"hair"); ii) the tangential arcs seen in the stellar distribution to 
the south ; iii) the ca 50 kpc long tail of atomic gas to the south. 

4.1. The model 

We have used numerical simulations to follow the merging pro- 
cess of two galaxies, a gas-rich disk-like galaxy (of spiral type) 
falling into a larger spherical system (an "elliptical" galaxy). 

The spherical galaxy is modeled by a set of collisionless 
particles distributed according to the Hernquist profile ( |1990| , 
with a potential of the form (]>(r) = -GMjir + a), where G is 
the gravitational constant, M is the mass, r is the radius and 
fl is a scale length. That model has a simple density profile 
(p oc r"'(r + aY^) and produces a surface density profile in 
good agreement with the de Vaucouleurs 7?'^"' law observed in 
elliptical galaxies, where R is the projected radius on the sky. 
Particles are assigned random velocities calculated from the en- 

{T990|. 



ergy distribution function given by Hernquist 



For simplicity, the elliptical galaxy is represented by a sin- 
gle component (collisionless particles in a single Hernquist 
model) and no distinction is made between stars and dark mat- 
ter Observationally, it is difficult to trace the mass distribu- 
tion in elliptical galaxies out to large radii and although the 
amount of data is rapidly increasing (e.g., Napolitano et al. 



2005 1 the scatter is large and the fraction of dark matter in 



elUpticals is stiU subject to interpretation (e.g., Romanowsky 



[eraL][2003] |Dekel et al.][2005 ^. No measurements of the stel- 
lar kinematics that would constrain on the mass-to-light ratio 
exist for the Medusa. In our model, the total mass of the ellip- 
tical is four times that of the infalling spiral. It corresponds to 
6.66 X 10"Mq, and about 50% of that mass is contained within 
a radius of 12 kpc. In the following we call the collisionless par- 
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tides "stars", but it should be kept in mind that some of them, 
especially at large radii, could be regarded as dark matter par- 
ticles. Since both the stars and dark matter particles have the 
same dynamics, having a mix of stars and dark matter woudn't 
affect the orbit of the infallling companion or the evolution of 
the system. 

The companion galaxy consists of a stellar disk, a gaseous 
disk, and a dark matter halo. The particles in the disk 
are distributed in two potentials of the form 0(R,z) = 
-GM/ -^^^2T|r(^^~^~y^r^^I)2^ where M is the mass, R is the 
radius in the plane of the disk, z is the coordinate perpendicular 



to the plane, and a and b are scale lengths ( Miyamoto & Nagai 



|1975 ). The advantage of that potential is that the scale lengths 
can be adjusted to represent different types of disks with a more 
or less prominent central bulge. A Hernquist model is used to 
represent the coUisionless halo. Initially, the particles are as- 
signed circular velocities in centrifugal equilibrium with the 
gravitational potential, and velocity dispersions are given by 
Toomre's stability criterion ( Toomre|1964 l. 

Initial conditions were generated from the distributions de- 
scribed above using Matlab. The evolution of the systems was 
computed by the programme COCKTAIL (Collisions of Clouds 
Keep Tree Algorithm in Labor). The coUisionless part is the 



tree algorithm of Barnes & Hut (1986i. The gas is modeled 
as an ensemble of clouds, all of the same mass, that dissipate 
energy through inelastic collisions with each other ("sticky par- 
ticles", e.g.. Combes & Gerin|1985||Kojima & Noguchi|1997| l. 
The gas accounts for about 10% of the mass of the stellar disk 
in the infalling galaxy. 

We used a total of 120000 particles: 30000 particles in the 
elliptical galaxy and 90000 particles in the spiral (40000 in the 
stellar disk, 20000 in the gaseous disk, 30000 in the dark matter 
halo). 

The parameters of the simulation presented here are listed 
in Table Setting the units of length and time to 1 kpc and 10^ 
years, with the gravitational constant G set to 1, gives a unit of 
mass of 2.22 x 10"^ Mq. We considered encounters of galaxies 
of different mass ratios and finally selected a mass ratio of 4:1, 
where the more massive galaxy is the elliptical. The initial mass 
of the gas in the spiral galaxy was set to match the observations 
of the Medusa, about 4.4 x IO^Mq. The mass of a gas particle 
is therefore 2.2 x 10^ Mq, which is of the order of that a giant 
molecular cloud. The mass of a star particle in the disk galaxy is 
10^ Mq, and it is about 22 x 10^ in the elliptical galaxy. Ideally, 
all particles of a given type should have a comparable mass to 
prevent artificial A^-body scattering (Chandrasekhar|[T942 1. In 
practice, however, this requirement is difficult to meet because 
of the limitation on the total number of particles set by an af- 
fordable computing time. 

Being mainly interested in the kinematical evolution of 
the stars and the gas in the infalling galaxy (and not in the 
internal relaxation of the merger remnant, which would be 
more sensitive to this issue), we used a comparatively larger 
number of particles to represent the stellar and the gaseous 
components of the infalling galaxy. 



Mass ratio 4: 1 

Elliptical: mass and scale length Meu = 300, a^n = 5 
initial position and velocity x = +10,y = 0,z = 

vx = 0, vy = 0, vz = 

Disk galaxy: Msiellardisk = 18, aslcllardisk = 1 

mass and scale length /'steiiai disk = 0.1 

-^gaseous disk 2, ^gaseous disk ^ 
^gaseous disk 0.1 
-'t^halo =55, flhalo = 5 

initial position and velocity x = -40, y = 0, z = 

vx = 0, vy = +0.2, vz = 

Table 4. Parameters of the simulations, in computing units. The 
units can be translated in the following physical units: 1 kpc, 
10^ yr, 100 km s"', 2.22 x IO'^Mq. 



The gas was initially distributed in a Miyamoto-Nagai disk 
with a scale length twice as large as that of the stellar disk, and 
the gaseous disk was truncated at an outer radius of 10 kpc. 
The scale length of the elliptical galaxy was calculated from the 
measured extent of the Medusa at the 25th magnitude/arcsec^ 
isophote: 025X1^25 = 2;3xl!6' - 26x18 kpc at our adopted dis- 
tance (RC3). This gives a radius to the 25th isophote of about 
1 1 kpc. The scale length of the Hernquist profile is related to 



the effective radius by the relation Rg = 1.8153a (Hernquist 
1990[). For elhptical galaxies, R^ < R25 < IR^. We thus took a 



scale length of 5 kpc for the elliptical galaxy. 

The two galaxies were initially separated by 50 kpc, and 
we tested different initial velocities. We obtained a best match 
for a slightly prograde encounter. In the simulation presented 
here, the spiral has initially a small positive tangential velocity 
(20kms-i). 

4.2. Results 

Figures[8]and|9]display the evolution of the stars and the gas as 
the spiral galaxy falls into the larger elliptical. Tidal distortions 
appear after closest passage and develop into a tidal tail, made 
of both stars and gas. The tail lengthens as the core of the in- 
falling galaxy and the larger elliptical galaxy orbit around their 
center of mass, leading to the eventual merging. Stars and gas 
are spread in the region diametrally opposite to the tail, form- 
ing a diffuse distribution remiscent of the "hair" of the Medusa. 
As particles from the spiral oscillate, shells form, due to the ac- 
cumulation of particles at their apocenters. 

The central part of the infalling disk galaxy is torn during 
the interpenetrating encouter with the more massive spherical 
galaxy. The material that will form the "hair" crosses the cen- 
ter of the spherical galaxy twice and starts expanding into the 
diffuse "hair" to the North while on the other side material ex- 
pands into the tidal tail. The core of the infalling galaxy falls 
back through the center of the spherical a third time, and that 
is when the simulation resembles the observations most. Stars 
from the companion oscillate back and forth and shells and 
loops are seen, but they are transient and their evolution can- 
not be followed in detail because of the relatively low number 
of particles that they contain. 
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Figure[TO]shows the distribution of the stars and the gas and 
the kinematics of the gas at the time when the best match with 
the observations was found, T - 355 Myr after the beginning 
of the simulation, after a rotation of +55° in the plane of the 
galaxy. At that time, the length of the tail is about 50 kpc, and 
the kinematics are in reasonable agreement with the observa- 
tions. Figure [To] can be compared both qualitatively and quan- 
titatively to Fig. |4] that shows the observed counterparts. The 
observed Hi velocity field shows no gradient along the tail, as 
in the simulation. The velocity dispersion map is also shown. 
The highest values are found in the central body, whereas the 
tidal tail shows little velocity dispersion. The isocontours show 
levels and constant particle density per pixel, where the pixel 
size of the simulated map is (1.5 kpc)^. A density of 100 parti- 
cles/pixel corresponds to 10^ Mq kpc"^, or a column density of 
hydrogen of 1.25 x 10^' cm"^, which translates into an Hi sur- 
face brightness of 0.5 Jy km s"' beam"'. In the simulation, the 
levels of highest gas density are therefore about twice as high 
as the highest contour in the observed map. This was to be ex- 
pected, since our total gas mass is exactly twice as large at the 
Hi mass and includes the molecular component, which is lo- 
cated in the central region of the galaxy. The simulation doesn't 
treat differently the atomic and the molecular gas. In the tail, 
the levels are comparable in the simulation and in the observed 
Hi map. Defining the tail as the region with y < -10 kpc (see 
Fig.[TO)i, about 30% of the gas particles and 20% of the star par- 
ticles lie in the tail. The gas fraction in the tail is comparable to 
the observed one. The amount of stars correspond to a surface 
brightness below the sensitivity of the optical observations. 

Let us now examine the kinematics of the gas in the simula- 
tion. The particles forming the extended tidal tail to the South 
and those in the diffuse "hair" to the North can be identified 
in the simulation by their different kinematics, as illustrated in 



Figure 1 1 The left figure shows the radial oscillations of the 
gas particles in the infalling galaxy. The orbit crowding result- 
ing in the tidal tail is clearly visible, starting at T = 210 Myr. 
The pattern is reminiscent to that due to radial waves in col- 
lisional ring galaxies (e.g. [Appleton & Struck-Marcell||I996| 
Horellou & Combes|200I i. Here, however, the phenomenon is 
different because there is no quasi-perpendicular impact with 
respect to the plane of the disk galaxy that would cause an ex- 
panding ring wave. The crowded particles seen in Fig. 1 1 are 
not located at all azimuthal angles, but come from an unwind- 
ing arm that forms in the infalling galaxy after the crossing of 
the elliptical. At T - 355 Myr, when the system resembles 
the Medusa most (the straight vertical line in the figure), most 
particles at a radius larger than 30 kpc are on their way out, 
wherease the ones at a smaller radius are falling back toward 
the galaxy. This is seen in a different way in the figure to the 
right, which displays the radial velocity of the gas particles ver- 
sus their radius. This is related to the slope of the curves in the 
left figure. Particles can be separated into three groups: the tail, 
expanding beyond r ^ 30 kpc {vrad > 0), the "hair", expand- 
ing beyond r ^ 12 kpc, and the rest, mostly located within 
r < 10 kpc. 

Since the different components have clear different kine- 
matics, we would like to be able to identify them in the ob- 
served Hi data. Unfortunately, since the system is seen very 



Table 5. The measured continuum flux, star formation rates 
based on continuum and FIR emission, respectively, and q pa- 
rameter 



Properties 


20 cm (Jy) 


0.108 


SFR20cm (Moyr-') 


22.8 


SFRFiR(Moyr-') 


9.5 


q 


2.5 



close to face-on, the kinematical signatures which appear in the 
radial plots are not present when we examine the line-of-sight 



velocity components, which are displayed in Figure 12 for both 
the simulation at T = 355 Myr (left) and the observations. The 
two plots are similar, displaying the gas tail extending to about 
60 kpc with a small and constant velocity dispersion as a func- 
tion of radius. The simulation shows more gas in the center 
than the observation which is sensitive to the atomic gas only, 
whereas the initial amount of gas in the simulation was taken to 
match both the molecular component (located in the inner part 
of the galaxy and traced by CO observations) and the atomic 
component. 

5. Discussion 

5.1. Star Formation Rate 



We foUowed |Condon| ([1992} and [Haarsma et al] pOOOj ) and 
used the relation below to calculate the SFR from the 20 cm 
continuum flux: 

SFR(Moyr-i) = 0.14D2F20cm 

where D is the distance in Mpc and Faocm is the 20-cm radio 
continuum flux density in Jy. The results are given in Table[5] 
We derived a SFR which is comparable to the results of |Condon 
et al. ( 1990| , who measured a flux of 122 mJy with the VLA 
and thus determined the SFR in the Medusa as ~ 20 Mq yr 

Table [5] also gives the SFR estimated from the FIR lumi- 
nosity: 

SFR(Moyr-i) = 0.17LFiR 



Kennicutt 



(1998 1, with Lfir in units of IO-^Lq. As 



following 

a reference for the IRAS 60/vm and 100/im flux densities in 



TabJTlwe used the IRAS Faint Source Catalog ( [Moshir &~et 
[aL] ( |1990} ). The relation between the SFR estimated from the 
20 cm flux and the FIR fluxes can be expressed by the q pa- 

(|1985 1, which is the logarith 



rameter defined by Helou et al. 



mic ratio of FIR to radio flux. Helou et al.| ( |1985| l derived a 
mean value of 2.3 for disk galaxies, independent of their star 
formation activity. Although NGC4194 is a merging galaxy, it 
conforms to this value which is a good indicator that the con- 
tinuum emission is dominated by emission from the starburst. 



Helou et al. ( 1985 i mention that even galaxies with a strong 
nuclear starburst, as Arp 220, are in agreement with this corre- 
lation of the FIR and radio flux, however within a higher scatter. 
Arp 220, e.g., has a q value of 2.65, whereas other prominent 
ULIRGs Uke Arp 299 and Mrk231 have a value of 2.1. 
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Our results are in good agreement with the investigations of 



Weistrop et al. (2004 1, who used the UV flux of the central star 
forming knots to estimate the star formation rate. They estimate 
a SFR within the knots of 5 - 6Mo yr ' and suggest an overall 
SFRofupto SOMoyr-i. 



The continuum source was studied recently by Beswick 
|et al.| ( |2005| l with sub-arcsecond resolution MERLIN data. 
They found a pair of continuum sources separated ~ 0.35", 
thus unresolved in our image. They derived a SFR of 6 Mg yr" ' , 



in agreement with the results of Weistrop et al. (2004|l, but they 



also pointed out that they only measured the star formation rate 
in the very center of NGC 4194. 

Our measured continuum source size of 1.2 kpc is only 
slightly smaller than the star forming region of ~ 2 kpc, as de- 
termined by optical and CO observations ( |Armus et al.||1990] 
lAalto & Huttemeister|20"00l l. 



5.2. The merger remnant 

Our simulations strongly support the hypothesis that the 
Medusa is the result of a merger between a gas-rich disk galaxy 
and a larger elliptical, with a mass ratio of about 4. We have 
tested a wide range of parameters, and the constraints placed 
by the observations of the Hi tail (its length and kinematics) 
and the optical appearance of the Medusa made it possible to 
reduce considerably the parameter space. The encounter that 
we finally selected is slightly prograde, in the sense that the or- 
bital plane coincides with the plane of the infalling disk galaxy, 
but the tangential velocity component is low - 20km s"'- and 
the infalling galaxy acquires largely radial motion as it falls 
through the larger elliptical. The two galaxies interpenetrate 
and the central part of the disk galaxy passes through the cen- 
ter of the elliptical twice, forming the diff'use "hair" on one 



side. Our simulations are similar to those of Kojima & Noguchi 
\991\ who investigated the formation of shells using a model 



similar to ours, where the gas is represented by "sticky" parti- 
cles. (Their model, however, included also star formation based 
on the local gas density and gas consumption). They showed 
that radial and slightly retrograde mergers lead to the formation 
of long-lived (> 1 Gyr) shells, and that star formation is sig- 
nificantly reduced because the gas is scattered in tidal debris. 
We have observed similar effects in some of the simulations 
that we have run. We believe that the Medusa is due to a dif- 
ferent type of encounter, namely a slightly prograde one, also 



discussed by 'Kojima & Noguchi ( 1997 1, that leads to both the 
formation of an extended gas-rich tidal tail and to the infall of 
gas toward the center, where a moderate starburst is observed. 
Although the overall distribution of the stars and the gas in the 
simulation is similar, some differences are seen. In low-density 
regions, the gas behaves like the collisionless stellar compo- 
nent. During the interpenetrating encounter, the rate of colli- 
sions between the gas clouds increase, which lead to a loss of 
energy and angular momentum and an infall of gas toward the 
center; the gas and the stars have a different initial radial distri- 
bution, with the gas being located more in the outer part with 
a larger angular momentum and lower binding energy. The gas 



is therefore more susceptible than the stars to be expelled into 
extended tidal tails. 

IHibbard et al. (j2000l studied some examples of disk- 
disk mergers which show anticorrelations between stellar and 
gaseous tidal features. They discuss various possibilities for 
these differences, including ram pressure stripping, photoion- 
isation, dust obscuration and different origins in the progenitor 
galaxies. Although possibly at least several scenarios may play 
a role at some point, the lack of starlight associated with the 
large Hi tail found in NGC 4194 can be easiest explained by 
the original position of the Hi in the progenitor 

6. Conclusions 

We have presented new optical, radio continuum and Hi line 
observations of the Medusa galaxy and numerical N-body sim- 
ulations of the merger of a small gas-rich disk galaxy with 
a larger elliptical galaxy. The optical observations consist of 
a deep V-band image of the distribution of the stars in the 
Medusa, and the radio observations of maps of the 20 cm ra- 
dio continuum emission and of the large-scale distribution and 
kinematics of the atomic hydrogen gas. The simulations are 
completely self-consistent with a total of 120000 particles. The 
gas in the infalling disk galaxy is represented by an ensemble 
of clouds, all of the same mass, that can dissipate energy in 
inelastic collisions with each other. 

1. The Hi observations revealed the presence of a single Hi 
tail to the South, extending out to a projected distance of 
4.95 ', or 56 kpc from the center of the Medusa. The tail 
contains ~ 7 ■ 10** M© of Hi and has no detected optical 
counterpart. No significant velocity gradient is observed 
along the tail, indicating that it is seen almost face-on. 

2. Hi was detected in a dwarf galaxy identified at a projected 
distance of 91 kpc to the North-West of NGC 4194. Its sys- 
temic velocity of2380kms"' is only 120 kms"' lower than 
that of the Medusa. This dwarf galaxy harbours ~ 7-10^ M© 
of Hi . No connection, neither in the optical nor the Hi 
could be found between the dwarf galaxy and NGC 4194 
and it is unlikely that it has had an influence on the dynam- 
ical evolution of the Medusa. 

3. We found a large misalignment between optical and 
gaseous tidal features. The diffuse optical tail is situated 
on the opposite side of the Hi tail, the optical shells do not 
seem to be connected to the Hi tail. Although the shells are 
within the region of the Hi tail, no connection is apparent 
neither in morphology nor in velocity distribution. 

4. We estimated a SFR of 23 M© yr"' based on the 20 cm con- 
tinuum flux density, and of 10 Mq yr"' from the FIR fluxes. 
This is consistent with a rather intense ongoing burst of star 
formation in the Medusa, in agreement with the optical and 



UV observations of Weistrop et al. (2004i 



Hi was detected in absorption toward the radio continuum 
in the center. 

The simulations of a slightly prograde merger between a 
gas-rich disk galaxy with a total mass four times less than 
that of the elliptical give the best match to the observations. 
A extended single tidal tail develops in the outer part while 
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the small galaxy falls in and oscillates around the center 
of the large elliptical, creating shell-like structures around 
the elhptical, best visible in the stars. The overall distribu- 
tion of the gas is similar to that of the stars, although more 
gas falls toward the center because of the dissipative na- 
ture and more gas is found of the extended tail because of 
the initially more extended distribution of the gas. The tail 
is mostly formed by material of the outer Hi disk of the 
progenitor spiral galaxy, and the siu'face brightness of the 
stellar material in the tail is low. 
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Fig. 4. Maps of the Hi distribution, velocity field and velocity dispersion, (top left) Hi distribution overlayed on optical image 
from the Digitized Sky Survey. The contour levels are 0.01,0.05,0.01,0.15,0.2,0.3,0.4,0.5 Jy beam km s^' (0.01 Jy/beam km/s 
corresponds to 2.47 • 10'^ cm"^), . (top right) Hi distribution as map and contours. Levels are the same. The white dot marks the 
region where Hi absorption is foimd. (bottom left) Velocity field, contour levels are from 2460 km s"' to 2600 km s"' in steps 
of 20km s (bottom right) Hi velocity dispersion, contour levels are 5, 10, 15, 20, 30, 40, 50km s"'. The synthesised beam is 
marked in the bottom left corner. 
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Fig. 8. Evolution of the stars in the simulation. The stars in the elliptical are shown in red, those in the infalling disk galaxy in 
black. 
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Fig. 10. View at a time corresponding to T = 355 Myr after the beginning of the simulation, when the system resembles the 
Medusa most. The orbital plane (which is the same as the plane of the disk of the infalling galaxy) also coincides with the plane 
of the sky. The system was rotated by 55° compared to the previous figures, a) Distribution of the stars, b) Distribution of the gas 
density. The grey scale shows the number of particles per 1.5 kpc^ pixel, c) Velocity field of the gas. The grey scale shows the 
range of velocities in km s~'. d) Map of the velocity dispersion of the gas. The grey scale shows the range of velocity dispersions 
in km s 
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Fig. 11. Kinematics of the gas in the simulation. Left: Temporal evolution of the radial distribution of the gas particles in the 
infalling galaxy. The vertical line at T = 355 Myr indicates the time when the system resembles the Medusa most. Right: Plot at 
T - 355 Myr of the radial velocity of the gas particles versus their distance from the center of mass of the infalUng galaxy. The 
"hair" and the tail appear as distinct kinematical features. 
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Fig. 12. Line-of-sight velocity of the gas versus distance from a center in the simulation (left) and in the observations (right). 



